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The response of cone photoreceptors to light is stable and reproducible because of the exceptional regulation of 
the cascade of enzymatic reactions that link visual pigment (VP) excitation to the gating of cyclic GMP (cGMP)- 
gated ion channels (cyclic nucleotide-gated [CNG]) in the outer segment plasma membrane. Regulation is 
achieved in part through negative feedback control of some of these reactions by cytoplasmic free Ca 2+ . As part of 
the control process, Ca 2+ regulates the phosphorylation of excited VP, the activity of guanylate cyclase, and the 
ligand sensitivity of the CNG ion channels. We measured photocurrents elicited by stimuli in the form of flashes, 
steps, and flashes superimposed on steps in voltage-clamped single bass cones isolated from striped bass retina. We 
also developed a computational model that comprises all the known molecular events of cone phototransduction, 
including all Ca-dependent controls. Constrained by available experimental data in bass cones and cone transduc- 
tion biochemistry, we achieved an excellent match between experimental photocurrents and those simulated by 
the model. We used the model to explore the physiological role of CNG ion channel modulation. Control of CNG 
channel activity by both cGMP and Ca 2+ causes the time course of the light-dependent currents to be faster than 
if only cGMP controlled their activity. Channel modulation also plays a critical role in the regulation of the light 
sensitivity and light adaptation of the cone photoresponse. In the absence of ion channel modulation, cone 
photocurrents would be unstable, oscillating during and at the offset of light stimuli. 



INTRODUCTION 

Cone photoreceptors in the vertebrate retina respond 
to light with high sensitivity over a large dynamic range. 
They respond to changes in luminance between dark- 
ness and the maximum radiance measured on earth's 
surface under solar illumination, ^1.6 x 10 9 cd/m 2 . At 
signal threshold, thoroughly dark-adapted cones re- 
spond to light flashes that excite as few as 4—12 visual 
pigment (VP) (cone opsin) molecules per cone, or a 
continuous stream of 30-60 excited VP molecules per 
second. Yet, they also adapt and respond to small per- 
centage changes in intensity, even when background is 
a steady stream as large as 10 10 excited VP/s (Burkhardt, 
1994; Paupoo et al., 2000). Over the first six orders 
of magnitude above threshold, cones respond with 
constant contrast: flashes of a given intensity, mea- 
sured as a percentage of the background intensity, 
generate the same amplitude response regardless of the 
absolute magnitude of the background luminance 
(Burkhardt, 1994; Normann and Werblin, 1974; Normann 
and Perlman, 1979). 

The sensitivity, speed, and adaptation of the cone 
electrical response are stable and reproducible because 
of the exceptional regulation of the cascade of enzy- 
matic reactions that couple the absorption of photons 
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to changes in membrane current. Regulation arises 
both from the interaction of the enzymes of phototrans- 
duction with regulatory proteins and the effects of cyto- 
plasmic Ca 2+ . In cones, Ca 2+ controls the enzymatic 
activity of cone VP kinase (VPK) (Kawamuraetal., 1996; 
Wada et al., 2006; Arinobu et al., 2010) and of guanylate 
cyclase (GC) (Lolley and Racz, 1982; Koch and Stryer, 
1988; Duda et al., 1996; Baehr et al., 2007; Takemoto 
et al., 2009), as well as the cGMP sensitivity of the CNG 
ion channels (Rebrik and Korenbrot, 2004) . 

The regulatory function of Ca 2+ in phototransduction 
is a feedback control mechanism because the output of 
phototransduction, a decrease in cytoplasmic Ca 2+ , in- 
fluences the events that lead to the change in Ca 2+ in 
the first place. It is a negative feedback because the out- 
put is fed back in such a way as to partially oppose the 
input. Signal stability is a serious challenge in the de- 
sign of engineering systems that are controlled through 
negative feedback (Storey, 2004). Moriondo andRispoli 
(2003) considered the instability (oscillations) that could 
occur in rod photocurrents in the absence of proper 
design of the regulatory Ca 2+ negative feedback loop. 
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To better understand the highly regulated activity of 
the molecular events underlying transduction in cone 
photoreceptors, we develop and analyze a detailed com- 
putational model of this process. The pioneering math- 
ematical modeling of phototransduction by Tranchina 
and Sneyd in cones (Sneyd and Tranchina, 1989; 
Tranchina et al., 1991) and Forti et al. (1989) in rods 
has been refined and renewed over time to yield con- 
temporary computational models that, although not 
without controversy, have made our understanding of 
phototransduction more complete and exact (in rods: 
Pugh and Lamb, 1993; Hamer, 2000; Caruso et al., 2005; 
Hamer et al., 2005; Shen et al., 2010; in cones: Reingru- 
ber and Holcman, 2008; Soo et al., 2008). Here, we de- 
velop a model of cone phototransduction that explicidy 
includes Ca-dependent control of VP phosphorylation 
and CNG ion channel cGMP sensitivity, features that 
are not a part of previous models. The model is tested 
by fitting and predicting experimental photocurrents 
measured under voltage clamp in isolated photorecep- 
tors. It is then used to investigate the physiological role 
of the Ca-dependent modulation of CNG channels; we 
find that CNG channel modulation is critically impor- 
tant to attain the time course, sensitivity, and stability 
characteristic of cone phototransduction. 

MATERIALS AND METHODS 

Materials and retinal cell dissociation 

3-12-mo-old striped bass (Morone saxatilis) were received from 
The Center for Aquatic Biology and Aquaculture at UC Davis and 
maintained at 15 °C in a small aquaculture facility under a 14:10 
(L:D) light cycle. Animal upkeep and experimental protocols 
were approved by the UCSF Institutional Animal Care and Use 
Committee and met all requirements of the National Institutes 
of Health Office of Laboratory Animal Welfare and the Associa- 
tion for Assessment and Accreditation of Laboratory Animal 
Care International. 

Fish were dark adapted for 30-40 min and then sacrificed in 
darkness. All of the following manipulations were performed 
under infrared (IR) illumination, aided by IR-sensitive video cam- 
eras and monitors. An eye was enucleated and hemisected, and 
the retina was then gently separated from the eyecup under 
glucose-Ringer's solution. Square tissue pieces, ^5 mm to a side, 
were recovered from the dorsal retina and incubated at room 
temperature for 3 min in 5 ml of glucose-Ringer's solution con- 
taining 0.5 mg/ml of collagenase from Clostridium histolyticum 
and 1 mg/ml bovine hyaluronidase (Sigma-Aldrich), followed 
by 5 min in hyaluronidase alone. The tissue pieces were then 
rinsed extensively in pyruvate-Ringer's solution to remove en- 
zymes and exchange the glucose. They were suspended in 500 pi 
of pyruvate-Ringer's solution and mechanically dissociated by 
teasing and tearing with fine forceps. 

Voltage-clamped photocurrent recording in isolated bass 
single cones 

The cell suspension was immediately deposited on a concanavalin 
A-coated glass coverslip that formed the transparent bottom 
of a recording chamber. The chamber was held on the stage 
of an inverted microscope equipped with differential interference 
conuast-enhancement optics and operated under IR illumination 



with the aid of IR-sensitive video cameras and monitors. After 
10 min in darkness, the chamber was vigorously perfused with 
glucose-Ringer's solution. Single and twin cones lacking their 
nuclear and synaptic regions (Miller and Korenbrot, 199.3) re- 
mained firmly attached to the coverslip. The chamber was in- 
termittently perfused throughout an experimental session, but 
not at the time photocurrents were measured. 

Tight-seal electrodes were fabricated from aluminosilicate cap- 
illary glass (1.5 x 1.1 mm; 1724; Corning) and applied onto the 
side of the cone inner segment. After forming a giga-seal, whole 
cell mode was attained by sustained suction while holding mem- 
brane voltage at 0 mV. Membrane current drifted continuously 
toward an outward (positive) value. When this drifting ceased 
(5-15 s; +40 to +70 pA), holding voltage was shifted to —40 mV, 
where membrane current was near zero. This method yielded 
more stable recordings than attaining whole cell mode at —40 mV. 
Voltage-clamped membrane currents were measured with a 
patch-clamp amplifier (Axopatch ID; Molecular Devices). Ana- 
logue signals were low-pass filtered below 50 Hz with an eight- 
pole Bessel filter (Frequency Devices) and digitally acquired at 
1 KHz (Digidata 1322A and pClamp 9.2; Molecular Devices). All 
photocurrents reported and analyzed here were measured 
within 8 min of the moment whole cell mode was achieved. 

Solutions composition 

The Ringer's solution was composed of 143 mM NaCl, 5 mM 
NaHC0 3 , 1 mM NaHP0 4> 2.5 mM KC1, 1 mM CaCl 2 , 1 mM MgCl 2) 
10 mM glucose, 10 mM HEPES, and lx MEM amino acids and 
vitamins, pH 7.5 (osmotic pressure of 309 mOsM) . In pyruvate- 
Ringer's solution, glucose was replaced with 5 mM Na pyruvate. 

Tight-seal electrodes were filled with a solution composed 
of (in mM): 115 K gluconate, 20 K aspartate, 20 KC1, 4 MgCl 2 
(0.5 mM free), 1 GTPNa 3 , 3 ATPNa 2 , and 10 MOPS, pH 7.25 
(osmotic pressure of 304 mOsM). The solution with all compo- 
nents except MgCl 2 was first made free of multivalent cations by 
passing it over an ion exchange resin in its K + form (Chelex 100; 
BioRad Laboratories). MgCl 2 was then added. 

Photostimulation 

Two optical benches were assembled, each with narrow band 
interference (540 ± 5 nm) and neutral density filters to control 
light spectrum and intensity. Stimulation was timed with electro- 
mechanical shutters (Vincent Associates) controlled with a digital 
pulse generator slaved to the data acquisition system. Collimated 
light beams from each bench were combined in a cube beam 
splitter and then focused by a condenser lens assembly (0.4 
numerical aperture [NA] ) onto the end of a 2-mm diameter liq- 
uid light guide of similar NA. The other end of the light guide was 
positioned at the point source in an epi-illuminator attached 
to the inverted microscope. A circular aperture within the epi- 
illuminator was focused by the microscope objective (40x, oil, 1.3 
NA; UV40; Nikon) to obtain a 40-pm diameter circle in the plane 
of the cells in the recording chamber. The outer segment of a 
cone under investigation was placed in the center of this circle. 
Light intensity was measured at the position of the recording 
chamber with a calibrated photodiode (UDT Sensors). Photon 
flux per unit area was converted to excited VP molecules (VP*) 
using 1.91 pm 2 as the collecting area of the bass cone outer seg- 
ment illuminated on its side with 540-nm wavelength light (Miller 
and Korenbrot, 1993). 

Computational analysis 

The system of simultaneous ordinary differential equations used 
to describe the molecular events of phototransduction was solved 
numerically using the Runge-Kutta method with integration time 
of 1 ps (20-Sim software v 4.2; Controllab Products BV) . The software 
also allowed parameter optimization using the Newton-Raphson 



32 



Modeling photocurrents in cones 



search method to reduce the difference between simulated and 
experimental data using Euler integral mean square error mini- 
mization. Optimized fitting of specified functions to experimen- 
tal data were also executed using Origin software (version 8.1; 
OriginLab) . 

RESULTS 

Photocurrents activated by flash and step photo 
stimulation of bass single cones 

Bass single cones are large photoreceptors (outer seg- 
ments are ^15-um long, 7 pm in diameter at their base, 
and 3 um in diameter at their truncated distal end) that 
express a single visual pigment with peak absorbance 
at 542 nm (Miller and Korenbrot, 1993). Photocur- 
rents measured at room temperature under voltage 
clamp in a thoroughly dark-adapted single cone in 
response to 10-msec flashes or 2-s light steps are illustrated 
in Fig. 1. Flash-elicited photocurrents reached a peak 



and then recovered with a distinct undershoot at the 
end ("off response). Peak amplitude is a function of 
light intensity well described by an exponential satu- 
ration function, 

I{VP*) = max /^aA(l-exp (-VP */&)), ( L1 ) 

where max /pm* is the maximum value of the photocur- 
rent, VP* is light intensity expressed as the number 
of excited VP molecules in the cell, and k is an adjust- 
able parameter (cr = ln2& is the light intensity at which 
current amplitude is half its maximum value) . We de- 
fined a signal threshold as the intensity that elicits a re- 
sponse of peak amplitude equal to twice the root mean 
square of the dark current noise. Mean values of ™*Ipeak, k, 
and signal threshold measured in 15 cones are listed 
in Table 1 . 

Steps of light 2 s in duration elicited photocurrents 
that reached an initial peak and then sagged over the 
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Figure 1 . Photocurrents measured under voltage clamp at —40 mV in a dark-adapted bass single cone in response to 10-msec flashes or 
2-s steps of 540-nm light. (A) Responses elicited with flashes of intensities: 9, 18, 37, 91, 178, 372, 913, 1,933, 3,685, and 9,133 VP*. Peak 
amplitude increases with light intensity (right) in a manner well described by an exponential saturation function (Eq. 1.1; continuous 
line), with Ipeak = 18.8 pA and k = 260.1 VP*. The inset is a log-log plot to show that signal threshold was ~8 VP*. (B) Responses elicited 
with step stimuli of intensities: 65, 327, 651, 1,635, 3,266, 6,513, 16,369, 65,131, and 163,687 VP*/s. The current at its peak (•) is more 
sensitive to light than in the stationary state (V; right). Amplitude dependence on light intensity is well described by a Michaelis-Menten 
function (Eq. 1.2; continuous line), with Ipeak =21.8 pA and cr = 1,844 VP*/ s for the response at its peak, and Imax =21.8 pA and 
rj = 3,522 VP*/s for the response at 2 s. The inset is a log-log plot of the peak amplitude of the step response to show that signal thresh- 
old was 56 VP*/s. 
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next second or so to a steady-state value (Fig. 1). The 
steady state reached in the time scale shown, however, is 
only apparent; under continuing illumination, a much 
slower phase of current recovery follows in a time scale 
of tens of seconds (Perlman and Normann, 1998). At 
the termination of the light step, photocurrents rapidly 
recover and exhibit a brief "off' response. Photocur- 
rents generated by very high intensity steps (over 10 4 - 
fold above threshold) exhibit a distinct, slowly decaying 
rebound in the dark current after the light step ends 
(Fig. 1; response to 35,000 VP*/s and above). The peak 
and steady-state amplitudes of the step response are a 
function of light intensity well described by the Michaelis- 
Menten function, 



I(VP*) = T 



VP* 
VP*+a' 



(1.2) 



where Imax is the maximum value of the photocurrent, 
is the rate of excited VP molecules per second, and 
cr is the intensity at half-maximum response. Mean val- 
ues of Ipeak, a at the peak and in the steady state, and 
signal threshold for 1 1 cones are listed in Table 1 . The 
photosensitivity of the cone decreases as the photocur- 
rent sags from peak to steady state, a manifestation of 
light adaptation. The loss of photosensitivity is evident 
by the fact that the value of cr at the peak is about twice 
that at the stationary state (Table 1). 

A model of the molecular mechanism of light transduction 
in bass single cones 

We developed a computational model that successfully 
describes the outer segment dark current and its 
changes with flash and step illumination. The model 
is based on the currently known molecular events of 
the phototransduction cascade and is constrained, 
when possible, by experimental facts known specifically 
for the bass single cones or generalized from data for 
cones in other species. The model follows and refines con- 
cepts introduced in previous models of cone phototrans- 
duction (Sneyd and Tranchina, 1989; Tranchina et al., 
1991; Reingruber and Holcman, 2008; Soo et al., 2008) 



but is novel in two major respects: (1) it develops in 
detail the inactivation of VP* by phosphorylation and 
its regulation by arrestin and Ca-dependent visinin 
(Kawamura, 1993; Kawamura et al., 1996), and (2) it 
introduces Ca-dependent modulation of ligand sensitiv- 
ity in the CNG ion channels (Rebrik and Korenbrot, 
1998;Rebrik etal., 2000). 

Biophysical and biochemical processes in the dark 
cGMP metabolic flux. In the dark, the cytoplasmic con- 
centration of cGMP is constant because the rates of 
its synthesis by GC and hydrolysis by phosphodiester- 
ase (PDE) are the same. This unchanging nucleotide 
concentration determines the magnitude of the dark 
current. 

The catalytic activity of GC and PDE in the dark in 
bass single cone is experimentally known (Holcman and 
Korenbrot, 2005): 

'<" rk VPDE = darh VGC = 6.5 uM/s. (2 ' 1} 

This value was used in all simulations of dark current. 

PDE activity in the dark. The cone-specific PDE holoen- 
zyme consists of two catalytic a subunits and two inhibi- 
tory 7 subunits (PDE6; see Conti and Beavo, 2007) . PDE 
enzymatic activity follows conventional Michaelis-Men- 
ten kinetics: 



dark\j^ 



PDE 



day!: 



^cat 



PDE 



[cGMP] 



dark 



[cGMP] dark + cGMP K v 



(2.2) 



where ^Vpde is the PDE hydrolytic velocity, cGMP K m is the 
Michaelis-Menten constant for cGMP, kcat is the cata- 
lytic turnover rate per active PDE molecule, ^pde* is 
the number of active PDE molecules in the dark, and 
[cGMP] dark is the nucleotide concentration. To simulate 
dark current there is no need to independendy know 
kcat and U *PDE* because it is possible to determine ex- 
perimentally the value of the parameter § ia & , where 



PDE 



(2.3) 



TABLE 1 

Light sensitivity of dark-adapted bass single cones 
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Mean 


26.2 


9.9 


245.3 


22.6 


14.9 


614.2 


106.05 


SEM 


2.4 


0.8 


25.6 


1.5 


1.0 


74.6 


86.4 


Cell count (n) 


15 


15 


15 


10 


10 


10 


10 


Minimum 


16.5 


4.6 


101.9 
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11.1 


315 


607 


Maximum 


43.6 


16.6 


452.1 


30 
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1,100 


1,612 
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The value c K m is known for cone PDE (Table 2) 
(Gillespie and Beavo, 1988; Huang et al., 2004), as is 
darh VpDE (Holcman and Korenbrot, 2005). For each cell, 
the dark current is measured, and from this value the 
concentration of cGMP in the dark can be computed 
(see below, Eqs. 2.5 and 2.6). Knowing these experi- 
mental quantities allows for calculation of the value of 
fidart in each and every cone analyzed (Table 2). 



GC activity in the dark. GC is a membrane-bound enzyme 
(GGE and GGF; Pugh et al., 1997; Koch et al., 2002, 
2010) , whose activity is controlled by Ca 2+ concentration 
as mediated by the Ca 2+ -binding proteins GCAP1 and 
GCAP2 (Dizhoor et al., 1995; Gorczyca et al., 1995; 
Kachi et al., 1999; Koch, 2002; Baehr et al., 2007). The 
biochemical activity of GC in fish cones has been recentiy 
characterized (Takemoto et al., 2009), but we define its 



TABLE 2 





Values of model parameters that best simulated dark current in 
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Idark 
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2.5 f 
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hi J dark 
J Ca 
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uM/s 


14.2 ± 5.4 s 
0.34 h 


18 
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20.9 ± 12.9' 
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0.052 + 0.024 


18 
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11.5 ±3.9' 


18 


Dark Ca 2 * efflux 


OUt J dark 
J Ca 


Statistical 


uM/s 


14.2 ± 5.4 g 


18 




max 
J Ca 


Statistical 


p.V 


4.87±1.88 i 


18 




exc 
Kca 


Adjustable 


( .M 


0.019 + 0.009i 


18 



invariant parameters are taken from the literature, and their values are held constant in all cells and all simulations of dark and photocurrents. "Statistical" 
are parameters measured experimentally in each and every cone analyzed. "Adjustable" parameters are arrived at by optimized fit of simulated to 
experimental data and are constrained, when possible, by data in the literature. The values in this table are the mean ± SD. 
Experimental data {Holcman and Korenbrot, 2005). 
Experimental data (Huang et al., 2004). 

Experimental data calculated in each cell from its measured dark current (and hence dark cGMP concentration), PDE K m , and dark PDE activity. See 
Eqs. 2.2 and 2.3. 

Experimental data in carp cones (Takemoto et al., 2009), but calcium dependence as measured in rods (Koch and Stryer, 1988). 
'Experimental data (Picones and Korenbrot, 1992; Rebrik et al., 2000). 

g Computed from Eq. 2.8 since dark current is measured. Vcos is 0.19 pL, about one half of the geometrical outer segment volume. The outer segment is 
a truncated cone of approximate dimensions: base diameter, 7.5 um; tip diameter, 3.5 vim; length, 17 pm (Miller and Korenbrot, 1993). Ca 2+ influx and 
efflux are the same in the dark. 
''Experimental data (Ohyama et al., 2000). 

'Experimental values in aequorin-loaded rods are estimated to be C m = 24 uM and £i-10 (Lagnado et al., 1992). 

j The time constant of Ca 2+ clearance from the bass cone outer segment determined from simulations with these values is 40.0 ± 15.4 msec, comparable 
to the experimental value of 43 + 9.8 msec measured in the tiger salamander cone (Sampath et al., 1999). Time constant is the time it takes to reduce 
cytoplasmic free Ca 2+ from 0.4 to 0.147 pM (1/e). 
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Ca dependence from biochemical studies of rod outer 
segments because cone GCAP is lost from the isolated 
fish cone outer segments (Takemoto et al., 2009). Elec- 
trophysiological studies in bass single cones suggest that 
the Ca 2+ dependence of GC activity in these cells is simi- 
lar to that in rods (Miller and Korenbrot, 1994): 



dark 



Vgc 



V 



max 
GC 



dark r \ n OC 

Ca 



1 + 



Ca ) 



(2.4) 



where ***Vgc is the GC catalytic velocity in darkness, 
is the maximum catalytic activity, dark Ca is the cyto- 
plasmic free Ca 2+ concentration in the dark, GC K Ca is 
the Ca 2+ concentration at which the enzyme activity is 
half its maximum value, and n GC is a dimensionless 
parameter that denotes cooperativity. The values of Vgc* 
and GC K Ca are known experimentally (Table 2) (Koch 
and Stryer, 1988; Takemoto et al., 2009) . 



CNG ion channel activity. The amplitude of the outer seg- 
ment current is determined by the activity of cGMP- 
gated (CNG) ion channels and is given by: 



I(cGMP,Co) 



[cGMP] 



n CNG 



[cGMP] +K c GMp( Ca ) 



(2.5) 



where I(cGMP, Ca) is the outer segment membrane cur- 
rent, Imax is the maximum cGMP-dependent outer seg- 
ment current, ^2,500 pA in bass cones at —40 mV 
(Rebrik et al., 2000), [cGMP] is cGMP concentration, 
and n C NG is a dimensionless parameter that denotes co- 
operativity. In bass single cones, ncNG = 2.5 and the value 
of J^cGiMp' 0 "' is Ca dependent (Picones and Korenbrot, 
1992; Rebrik and Korenbrot, 1998; Rebrik et al., 2000): 



We assigned Ca = 0.4 pM based on experimental 
measurements in tiger salamander cones (Sampath 
et al., 1999). Under this assignment, we used Eqs. 2.5 
and 2.6 to compute free cGMP in each cone from its 
measured dark current (Table 2). 



Ca 2+ ionic flux. It has long been established that Ca 2+ ions 
flow into the photoreceptor outer segment via the CNG 
ion channels and flow out via Na + /Ca 2+ ,K + exchangers 
(Yau and Nakatani, 1985; Miller and Korenbrot, 1987) . 
In the dark, Ca 2+ influx and efflux are identical, and the 
outer segment free Ca 2+ concentration is constant. 



n J 



dark 
Ca 



t Td 
JC 



Ca 



(2.7) 



Dark Ca 2 * influx. Ca 2+ ion influx, m J a (pM/sec) , is a con- 
stant fraction of the inward outer segment current at all 
voltages (Ohyama et al., 2002) and is given by: 



in j 
JC, 



ZP/10 
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Ca 



zFV C0S Buff(Ca) 



(2.8) 



where I (pA) is the membrane current, Pf is the fraction 
of the current carried by Ca 2+ , % is the Ca 2+ valence, F is 
Faraday's constant, V cos is the outer segment cytoplasmic 
volume, and Bujj(Ca) is the cytoplasmic Ca 2+ -buffering 
capacity, defined below. The values Pf = 0.34 and V cos = 
0.19 pL for the bass single cone are known from experi- 
mental determinations (Miller and Korenbrot, 1993; 
Ohyama et al., 2000) . Dark Ca 2+ influx is computed from 
Eq. 2.8 with current, /, set at the value of the dark current. 

Dark Ca 2+ efflux. Ca 2+ ion efflux, mt J c a (pM/sec) depends 
on free cytoplasmic Ca 2+ and is given by: 
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where out J ^™ ax (in pA) is the maximum exchanger-me- 
diated efflux, and other terms are the same as in Eq. 2.8. 
The term in the square parenthesis is the Ca 2+ depen- 
dence of the Na + /Ca 2+ ,K + exchanger transport rate, 
where Ca is free cytoplasmic Ca 2+ in the dark and 
is the Ca 2+ concentration at which the transport rate is 
half its maximum value (Sheng et al., 2000) . 

Because we measured the dark current in every cell 
we studied, we could calculate the magnitude of Ca 2+ 
dark influx (Eq. 2.8) and, therefore, know Ca 2+ dark 
efflux (Eq. 2-9). Assigning Ca = 0.4 pM in darkness 
(Sampath et al., 1999) allowed us to adjust the values of 
On average in bass single cones, the values of mmK c GMP Jo, and ifc a (Tables 3-5). Although these two values 
and m! " K cGMP are 105.5 and 316 pM, respectively, and are indeterminate if dark current alone is simulated, re- 
CN ° K Ca = 0.86 pM (Rebrik and Korenbrot, 1998; Rebrik quiring that the same values simulate dark and photocur- 
et al., 2000) . rents in the same photoreceptor constrained their values. 



i mm T - i max T , • i 

where K C CMP an( * K cCMP are minimum and maxi- 
mum values of k c gmp - Ca is cytoplasmic free Ca 2+ con- 
centration, and CWG ^ca is the Ca 2+ concentration at which 
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Cytoplasmic Ca 2+ buffer. The molecular identity and 
physiological characteristics of the cone outer segment 
Ca 2+ buffer are unknown, but studies in rod outer seg- 
ments suggest that there exist two instantaneous buffer 
systems: one of high affinity and low capacity that oper- 
ates below 1 uM free Ca 2+ , and another of low affinity 
and high capacity that operates at higher Ca 2+ levels 
(Lagnado et al., 1992). Assuming that similar Ca 2+ 
buffers exist in the cone outer segment cytoplasm, 
and following Lagnado et al. (1992): 



Ca TotB = HA +(B + V)Ca, 
Ca + K HA 



(2.10) 



2+ 



where Ca TotB is the total exchangeable (bound) Ca 
Cha is the total capacity (concentration) of a high affin- 
ity buffer of Michaelis-Menten constant K m , and B is 
the buffer capacity of the low affinity buffer. 

Ca 2+ -buffering capacity, the ratio of bound to free Ca 2+ , 
in photoreceptor outer segments is not constant but 
changes with Ca 2+ concentration. To calculate the Ca de- 
pendent Ca 2+ buffer capacity, Buff(Ca), the derivative of 
Eq. 2.10, is computed (Berlin et al., 1994; Neher, 1995): 



argument that six is a reasonable number for rhodopsin, 
which we have adopted. Kennedy et al. (2004) have dem- 
onstrated that light causes phosphorylation of up to five 
different sites per cone VP molecule in vivo. 

The catalytic activity of VPK is regulated by its inter- 
action with two different proteins: (1) arrestin, which 
quenches VP* lifetime by binding phosphorylated VP* 
and precluding further G protein activation (Craft and 
Whitmore, 1995; Gurevich and Gurevich, 2006; Nikonov 
et al., 2008) ; and (2) visinin (recoverin in rods) , a Ca-bind- 
ing protein that mediates Ca-dependent regulation of 
VPK activity (Kawamura, 1993; Kawamura et al., 1996). 

VPK adds phosphates one at a time at a rate of 1n Pi > 
which depends on n m , the number of phosphorylated 
sites. Each phosphorylated state of VP* is identified as 
w n Pi > where o> n ft <6. When n Pi = 6, VP* is inactive. 
Arrestin binds to phosphorylated VP* with a rate ^n Pi 
and quenches the ability of VP* to activate transducin. 
Therefore, the rate of change in the number of VP* 
molecules with n Pj phosphorylated sites is given by: 



dW 0 (t) 
rtl 



Phot(t) - y 0 VP 0 (t) 



(3-1) 



Buff(Ca) -- 



dCa Tl 



otB 



dCa (Ca + K HA ) 



- + B + 1. (2.11) 



Simulations were very sensitive to the values of K m 
and Cat- To constrain these parameters, we assigned 
C m an initial value of 20 uM and B= 10, consistent with 
those measured in rods. We adjusted the value of to 
optimally fit simulated to experimental dark and photo- 
currents using a computer-aided least-square minimiza- 
tion search method (Raphson-Newton; tolerance, 0.001) . 
We demanded the same optimized values of these pa- 
rameters fit dark and light-dependent currents in the 
same cell. Mean values of these parameters arrived at by 
optimized fit of simulated to experimental membrane 
currents are listed in Tables 2-5. 

Biophysical and biochemical processes in the light 
Kinetics of VP activation and inactivation. Photon absorp- 
tion generates an active state of cone VP, VP*, one that 
interacts with the G protein transducin to initiate the 
phototransduction process (Hofmann et al., 2009). 

VP* inactivates because it is phosphorylated, a pro- 
cess catalyzed by the activity of a cone-specific kinase 
(G protein-coupled receptor kinase 7; VPK Hisatomi 
etal., 1998; Rinner etal., 2005; Tachibanaki etal., 2005). 
The time course of VP* inactivation is complex because 
VPK catalytic rate depends on the state of VP* phos- 
phorylation: the rate of phosphorylation of additional 
sites depends on the number of sites already phosphory- 
lated. The number of physiologically relevant VP* 
phosphorylation sites is controversial, as few as 2 and as 
many as 16. Hamer et al. (2005) have made a persuasive 



dVP l (t) 
dt 



^YoWoM-tYi+Hi)^ (t) 



(3-2) 



dVPnpit) 
dt 



n Pi 



(3.3) 



where Phot(t) is the number of VP* molecules produced 
by stimulus photons, Ywp,- is the rate of phosphoryla- 
tion of VP* with n pi phosphorylated sites, and M-n ft - is 
the rate of arrestin-dependent quenching of VP*. 

The rate of VP* phosphorylation ^n pi decreases expo- 
nentially as the number of phosphorylated sites increases: 



Y(n Pi +l) ^ n Pi' 



-CDy 



(3.4) 



where is a constant. ^n pi is assumed to increase lin- 
early with the number of phosphorylated sites n pj • 



^~ n Pi ~ npi ^ 0 ' 



(3.5) 



The Ca 2+ dependence of cone VP* phosphorylation rate 
was taken to be the same as that of the rhodopsin kinase 
(our fit to the descriptive data) (Kawamura, 1993): 



Yo ~ Jmax 



0.1 + 0.9 



1 + 



Ca 



IK, 



Ca J 



(3.6) 
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where Ymax is the maximum possible value of this rate, 
and y K Ca = 0.9 uM. 

cGMP metabolic flux 

Time course of light-activated PDE activity. VP* interacts 
with the G protein transducin and generates T*, an 
activated form of the G protein. The rate of T* creation 
is given by the encounter frequency between VP* and 
transducin in the plane of the membrane, a feature lim- 
ited by the diffusion of the colliding molecules (Lamb, 
1996) . ^nPi is a rate gain that specifies the number of 
T molecules activated per second by one VP*, a value 
that changes with the number of phosphorylated sites 
in VP* (Gibson et al., 2000). The affinity between of 
VP* and T declines exponentially with approximately 
twofold decrease per phosphorylation. Thus, 



(n Pi +l) 



-co. 



it 



act 



Pi 



(3.7) 



where, again, n Pi is the number of phosphorylated 
sites, between 0 and 6, and C0 ac ; = 0.69 is the rate of expo- 
nential decay, an experimentally known value (Gibson 
etal., 2000) (Table 3). 

T* activates PDE with a 1:1 stoichiometry (Leskov 
et al., 2000) to produce PDE*. PDE* inactivates at a rate 
a- pde- Therefore, 



dPDE * (t) 
dt 



l Pi 



n pi -UpDEPm*(t)- 



(3.8) 



T* disappears at a rate that is commonly thought to 
be the rate of GTP hydrolysis by GTPase activity inher- 
ent to Tct-GTP. This is not the case in rod photorecep- 
tors. Disappearance of T* requires the interaction of 
Ta-GTP with several other proteins: RGS9, R9AP, G(35L, 
and yPDE (Wensel, 2008). Recent analysis has shown 
that the kinetics of this multi-molecular interaction is 
well represented by a simple "Michaelis module" in 
which the multiprotein complex behaves as though a 
single protein, "RGS9," interacts with T*; the rate of for- 
mation of this T*-RGS9 complex is slower than the rate 
of GTP hydrolysis (Burns and Pugh, 2009) . In such a 
scheme, then, the aPDE term of Eq. 3.8 is the rate 
of interaction of the RGS9 complex with T*, not the 
rate of GTP hydrolysis. An important consequence of this 
scheme is that the value of oipde becomes dependent 
on RGS9 concentration, an experimental fact (Krispel 
etal., 2006). 

Light-dependent changes in cytoplasmic free Ca 2+ and the 
enhancement of GC activity. Cytoplasmic free Ca 2+ 
changes when its rate of influx is different than its 
rate of efflux (Yau and Nakatani, 1985; Miller and 
Korenbrot, 1987). Therefore, the light-dependent rate of 
change in Ca 2+ concentration is given by the difference 



of Eqs. 2.8 and 2.9, as membrane current, I(t), changes 
with light: 



(3.9) 



I(t)Pfl0 6 -2" ut J™ K 



dCa(t) 
dt 




zFV cos Buff(Ca) 



Integration of Eq. 3.9 yields the magnitude and time 
course of light-dependent changes in free Ca 2+ concentra- 
tion. Light-dependent changes in GC activity simply track 
the changes in Ca 2+ concentration, as defined in Eq. 2.4. 

Time course of light-dependent changes in cytoplasmic 
cGMP concentration. The light-dependent rate of cGMP 
hydrolysis is determined by PDE* catalytic activity, 
'• gi 'V , where: 



llgk % DE (VP*,cGMP) = $ sub PDE * -g^g 



(3.10) 



cGMP(t) + cGMF K v 



PDE* is the number of active PDE (the integral of 
Eq. 3.8), cGMP is the nucleotide concentration, cGMP K m 
is the PDE Michaelis-Menten constant for cGMP, and 



2^10" 



(3.11) 



p sub is used because PDE* is in units of number of mol- 
ecules, k cat is the enzymatic turnover rate per active 
PDE, and 2 reflects the fact that an active PDE* is a dimer. 
L is Avogadro's number, and Vcos is the cone outer seg- 
ment cytoplasmic volume. In simulations, the value of 
P sl/ j is adjusted. 

In the course of the photoresponse, GC activity 
changes with time because cytoplasmic Ca 2+ changes 
(the integral of Eq. 3.9): 



R Z h % c (Ca 



V 



I1KIX 
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1 + 



Ca(t) 

KCa 



(3.12) 



GC 



The net rate of change of cGMP in the course of a light- 
elicited response is: 



dcGMP(yP*,Ca,cGMP) (da * v 

= ( Vac +V GC (Ca))- 

at 

C ark VpDE + " g '"V PDE (VP*,cGMP)). 



(3.13) 



cGMP-gated, Ca-modulated membrane current. Integra- 
tion of Eq. 3.13 yields the time course of light-dependent 
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changes in cGMP. The simultaneous changes in Ca 2+ 
and cGMP concentration are then used to compute 
the light-dependent changes in membrane current 
(Eqs. 2.5 and 2.6). 

Optimized fit of model simulations to experimental data 
Selection of parametric values. To fit simulated to experi- 
mental data, we divided model parameters into three 
categories: invariant, statistical, and adjustable. "Invari- 
ant" parameters were selected from the available lit- 
erature, and their values were kept the same for all 
cells and all simulations. "Statistical" parameters were 
experimentally measured in the cell under analysis. 
"Adjustable" parameters were adjusted to optimize 
the fit between simulated and experimental currents in 
each cell under study. Adjustments were made with 
computer assistance to fit simulated to experimental 



data by least-square minimization (Raphson-Newton; 
tolerance, 0.001). 

The values of adjustable parameters were not arbi- 
trary; their initial values were constrained, whenever 
possible, by experimentally known values available for 
bass single cones specifically or from biochemical studies 
of cone phototransduction. Values of adjustable parame- 
ters are presented in Tables 2-5, and footnotes quote the 
constraining experimental values. Simulations were par- 
ticularly sensitive to the following adjustable parameters: 
(a) The features of the Na + / Ca 2+ ,K + exchanger transport, 
°"'Jca X and K^, constrained by demanding their value 
be the same in the dark and the light for the same cell. 
Furthermore, their values in the dark were constrained in 
each cell because the dark Ca 2+ efflux was known (from 
the dark current; Eqs 2.7 and 2.8), and cytoplasmic free 
Ca 2+ in the dark was set at 0.4 pM (Sampath et al., 1999; 



TABLE 3 

Values of model parameters that best simulated photocurrents in dark-adapted bass single cones: Response to 10-msec light flashes 





Parameters 


Category* 


Units 




Cone 1 




Cone 2 


Cone 3 




Mean ± SD b 


VP* 


Intensity 


Statistical 


VP* 


167 


17,443 


174 


17,136 


173 


7,720 


Fig. 4 




i max 


Adjustable 


uM/s 


100 c 


65 


100 


70 


105 


6S 






Vh) 


Invariant 


1/s 


0.5 


0.5 


0.5 


0.5 


0.5 


0.5 






°\ 


Invariant 




0.1 


0.1 


0.1 


0.1 


0.1 


0.1 








Adjustable 




230 d 


230 d 


215 


215 


230 


230 


231 + 19 




®acl 


Invariant 




0.69 e 


0.69 e 


0.69 e 


0.69 e 


0.69 e 


0.69 e 




PDE* 




Adjustable 


pM/s per 
molecule 


0.185 


0.185 


0.323 


0.323 


0.416 


0.416 


0.334 ± 0.091 f 




®-PDE 


Adjustable 


1/s 


17 B 


47 


7 


10 


14 


28 


Fig. 4 


CNG ion 


Mark 


Statistical 


pA 


22.2 


22.2 


18.8 


18.8 


42.3 


42.3 


27.3 + 10.5 


channel 






















Ca 2+ buffer 


K HA 


Adjustable 


|iM 


0.030 


0.030 


0.044 


0.044 


0.091 


0.091 


0.052 + 0.025 




('HA 


Adjustable 


|i.M 


21.4 


21.4 


4.6 


1,6 


66 


66 


20.9 ± 12.9 




B 


Adjustable 




10.6 


10.6 


7 


7 


15 


15 


11.5 ± 3.9 


Ca 2+ efflux 


/"max 
JCa 


Statistical 


pA 


4 


4 


3.23 


3.23 


7.63 


7.63 


4.87 + 1.88 h 




kZ 


Adjustable 


|i.M 


0.025 


0.025 


0.005 


0.005 


0.025 


0.025 


0.019 ± 0.009 h 



a Values of invariant parameters first used to compute dark current are not repeated here. They have the same values listed in Table 2. 
b Mean ± SD; n = 18 single cones. 

c Imax experimental value is 100-150 s -1 (Tachibanaki et al., 2005). 

d This value was initially constrained to 150 s -1 from experimental data (Leskov et al., 2000). 
Experimental data (Gibson et al., 2000). 

'Experimental values of PDE* k cat range from 2.2 to 8 x 10 3 s -1 (Dumke et al., 1994; D'Amours and Cote, 1999; Muradov et al., 2009). From the definition 
of $sub (Eq- 3.11) and the highest k cat value, the 0.334 value indicates the effective intracellular volume for cGMP is, on average, 0.08 pL, rather than the 
0.19 pL used throughout these calculations. This factor of approximately two can probably be accounted for by the fact that the volume, in fact, changes 
some from cell to cell and is not the rigid assigned value of 0.19 pL. There is experimental uncertainty because the size of each and every cone recorded 
from was not measured. Also, the actual concentration of cGMP may be a little smaller than the computed value because there is intracellular cGMP 
binding (Olson and Pugh, 1993), a fact that would be revealed as a smaller apparent volume. 

g This value was initially constrained to 10 s -1 at dim lights from the analysis of dark current noise in bass single cones (Holcman and Korenbrot, 2005). 
h With these values, the computed time constant of Ca 2+ clearance from bass cone outer segment is 40.0 ± 15.4 msec, similar to the experimental value of 
43 + 9.8 msec in the tiger salamander cone (Sampath et al, 1999). 
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Leung et al., 2007). (b) The features of cytoplasmic Ca 2+ 
buffers, and K HA , constrained by demanding their 
value be the same in the dark and the light, and their 
initial values consistent with data available for Ca 2+ buff- 
ering in rod outer segments (Lagnado et al., 1992). Fur- 
thermore, the values were constrained by the need to 
simulate the dark current measured in each cell with 
free Ca 2+ set at 0.4 uM (Sampath et al., 1999). (c) lmax> 
the maximum rate of Ca-dependent VP* phosphoryla- 
tion, with initial value 100 s from measurements of 
VPK activity in fish cone outer segments (Tachibanaki 
et al., 2005). (d) a PDE , the rate of PDE* inactivation 
with initial value 10 s _1 from measurements of bass 
cone outer segment dark current noise (Holcman and 
Korenbrot, 2005). 

Both VP* phosphorylation ( y max ) and PDE* inactiva- 
tion ( a mE ) contribute to the control PDE* lifetime, and 
there is no a priori mathematical argument to decide 
whether one is slower than the other and, therefore, the 
dominant rate constant that controls the duration of the 
photoresponse (Burns and Pugh, 2010). However, cer- 
tain experimental facts are known for cones: (a) o- PBE is 



5-15 s _1 at dim lights, inferred from dark current noise 
(Holcman and Korenbrot, 2005) , much slower than y m ax 
(100-150 s _1 ) measured biochemically (Tachibanaki 
etal., 2005). (b) y m ax is the dominant rate constant that 
controls photocurrent duration at high light levels that 
saturate cone photocurrent amplitude and above, as 
demonstrated by electrophysiological experiments 
(Matthews and Sampath, 2010). (c) At light levels 
that bleach 0.1% of the VP and above, cone photocur- 
rent duration is controlled by the rate of decay of the 
cone-opsin photoproduct Metall, shown by electrophysi- 
ological measurements combined with VP chromophore 
substitution experiments (Estevez et al., 2009). Because 
of this evidence, we adopted the postulate that at intensi- 
ties below photocurrent amplitude saturation, ymax Was 
constant and only the value of a mE was adjustable, 
whereas for intensities at and above amplitude satura- 
tion, a mE was constant and only y ma x was adjustable. 
This hypothesis on the dual control of PDE* lifetime was 
not only consistent with available information but also 
proved critical for successful photocurrent simulations, 
as shown below. 
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Figure 2. Experimental (gray traces, noisy) and simulated (black traces, noiseless) photocurrents measured at —40 mV in a dark- 
adapted bass single cone. Photocurrents were elicited by 10-msec light flashes of intensity: 36, 71, 167, 356, 710, 1,744, 3,561, 7,106, and 
17,443 VP*. The values of the parameters used to compute the simulated data are listed in Table 3 (Cone 1) . 
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Responses of dark-adapted cones to light flashes 

Photocurrents measured in a bass cone in response to 
10-msec light flashes of varying intensity and simulated 
currents optimally fit to the experimental data are illus- 
trated in Fig. 2. Values of the model parameters used in 
these computation are presented in Table 3 (Cone 1). 
Simulated photocurrents fit experimental data well, 
except they did not fit the small exponential compo- 
nent observed as photocurrent amplitude approached 
saturation. This component is generated by the electro- 
genic activity of the Na + /Ca' + ,K + exchanger (Nakatani 
and Yau, 1989; Perry and McNaughton, 1991), and the 
model does not specifically include simulation of this 
active transport system. To report the range of quality 
of fit between model and experimental data, Fig. 3 
illustrates experimental and simulated data for two ad- 
ditional cones. Table 3 lists parametric values corre- 
sponding to simulations in Fig. 3, as well as the mean ± 
SD of values determined in a total of 18 cones. 

The values of the parameters o- PDE and Tmflx were 
found to be light dependent. Their mean value as a 
function of intensity arrived at by successful simulations 
in 18 cones is shown in Fig. 4 B. At all intensities below 



amplitude saturation, ^2,000 VP* PDE inactivation, 
o-pde was the dominant rate constant in control of photo- 
current duration. Its value changes little over the first 
^500 VP* and then rises up to a maximum value, where 
it remains over all intensities tested. At intensities above 
amplitude saturation, the rate of VP* inactivation by 
phosphorylation y max becomes the dominant rate con- 
stant. As is illustrated in Fig. 3 (B and D) , simulations fit 
experimental data only if the value of y m ax was assigned 
a rate-limiting role and adjusted. If the value of l max 
was not adjusted, but kept at the value used to fit non- 
saturated photocurrents, simulations failed to fit exper- 
imental photoresponses (Fig. 3, B and D, red lines) . 

It is possible to adjust o- PDE to fit the experimental 
photocurrents over all intensities tested. However, this 
would yield a biphasic dependence of o. PDE on light, at 
first increasing with light intensity and then decreasing 
as the stimulus becomes brighter. Such biphasic behav- 
ior seems unreasonable and without a basis of fact. 
Thus, the hypothesis that PDE* lifetime is controlled by 
a dual mechanism, the rate of PDE inactivation at low 
intensities and VP* inactivation by phosphorylation 
at intensities above amplitude saturation, is consistent 
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Figure 3. Experimental (gray traces, noisy) and simulated (black traces, noiseless) photocurrents measured at —40 mV in two different 
dark-adapted bass single cones. Data in A and B were measured in the same photoreceptor. Flash intensities tested were 174 and 1,747 
VP* in A and 1 7,743 VP* in B. C and D were measured in the same cell. Flash intensities tested were 1 73 and 808 VP* in C and 7,720 VP* 
in D. The values of the parameters used to compute the simulated data are listed in Table 3. At intensities above amplitude saturation 
(B and D), simulations fit experimental data only if the VP* inactivation rate, lmax> was slower than for dim light responses (Table 3). 
The simulations illustrated in red were obtained when "imax was kept at the same value used to fit the dim light responses. 
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with experimental findings in other cones (Matthews 
and Sampath, 2010) and allows for successful simula- 
tion of experimental photocurrents in bass cones. 

The Na + / Ca 2+ ,K + exchanger transport rate determines 
the velocity at which Ca 2+ is cleared from the outer 
segment when all light-sensitive channels are suddenly 
closed (Miller and Korenbrot, 1987). Experiments show 
the Ca 2+ clearance time constant in tiger salamander 
cone outer segments is 43 ± 9.8 msec (±SD; n = 25) 
(Sampath etal., 1999). The Ca 2+ clearance time constant 
computed from simulations using optimally adjusted 
values of /™ x and K e £ (Tables 2-5) is 40 ± 15.5.5 msec 
(±SD; n = 18), in good agreement with experiments. 

Based on experimental findings in rods (Lagnado etal., 
1992), we assumed the existence of both high and low 
affinity Ca 2+ buffers. The features of this buffer system, its 
capacity and dissociation constant, were arrived at by 
computer-aided optimal fit of simulated to experimental 
data, constrained by the experimental values reported 
for rod outer segments (Tables 2-5) (Lagnado et al., 
1992). Ca 2+ buffer capacity is a function of Ca 2+ (Eqs. 2.10 
and 2.11), and because Ca 2+ changes with light, buffer 
capacity is also light dependent. To illustrate this depen- 
dency, we measured the simulated free Ca 2+ at the peak 
of the photocurrent at each intensity and calculated buf- 
fer capacity using the mean values of B, c m , and ^ m in 
Table 2. Buffer capacity, dC ^ B > increases with light inten- 
sity (Fig. 4) , and its value is consistent with measurements 



of this quantity in many other cells, where values range 
between 40 and 480 (Neher, 1995). 

Responses of dark-adapted cones to steps of light 
In rods, photocurrents elicited by steps of dim lights are 
the temporal summation of individual dim flash re- 
sponses (Baylor, 1987). Similarly, cone photocurrents 
elicited by light steps under ^600-800 VP*/s (near the 
cr value of the initial peak of the step response) are well 
simulated by the temporal sum of dim flash photocur- 
rents. Fig. 5 shows simulated and experimental voltage- 
clamped photocurrents measured in one dark-adapted 
single cone in response to 2-s light steps of varying 
intensity. Responses to stimuli under 1,177 VP*/s were 
well fit by simulations with the same parametric values 
(Table 4, Cone 1 ) inferred from the simulation of flash 
photocurrents in the very same cone. The sum of dim 
light flash photocurrents did not simulate the responses 
to light steps brighter than 1,177 VP*/s. 

The need to evoke additional molecular processes to 
explain the sag in the step response from a peak to a 
stationary value was first recognized by Forti et al. (1989) 
in rods. They proposed that the sag arises from time- 
dependent changes in the concentration of a Ca 2+ -binding 
buffering molecule. This mechanism is improbable be- 
cause it predicts light-dependent changes in cytoplasmic 
Ca 2+ that are far different from what is observed experi- 
mentally. Nikonov et al. (2000) demonstrated that the 
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Figure 4. Light dependence of the values of adjust- 
able parameters inferred from successful simula- 
tions of flash photocurrents. Data points are the 
mean ± SD of values computed in 18 single cones. 
The light intensity shown on the graph is the mean 
of intensities binned within about plus or minus 
10% of the mean. (A) PDE* inactivation rate ( Ct PDE ) 
and VP* phosphorylation rate ( ''-max ) . (B) Ca 2+ - 
buffering capacity ( dCa " T ° um / d Q \ Eq. 2.11). Continu- 
ous lines are drawn by eye to join the data points. 
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sag in rods likely arises from a time-dependent increase 
in the PDE* inactivation rate. In cones, on the other 
hand, Soo et al. (2008) did not propose an a priori bio- 
chemical hypothesis to explain the sag in the step re- 
sponse; rather, they used simulations and found an 
optimal fit of simulated to experimental step responses 
by allowing an increase in PDE* inactivation rate (just as 
in rods) , as well as changes in the Ca dependence of GC 
activity and the Na + / Ca 2+ ,K + exchanger transport rate. 

In our simulation, in agreement with Nikonov et al. 
(2000) , we assumed that when summation of dim light re- 
sponses fails to fit experimental step photocurrents and 
up to intensities that saturate current amplitude, there is a 
time-dependent increase in PDE inactivation rate. Simula- 
tions successfully fit experimental photocurrents (Fig. 5) 
under the assumption that cl pde accelerates from an initial 
value m "a PDJ£ to a final value f '" al o. mE =""' a PBE + a<x pbe along 
an exponential time course. Thus: 



<W0= rt a m + Aa m (i- 



- h *PDE t )) (3-14) 
where 

pde i s tri e exponential rate of change. 



To report the range of quality of fit between model 
and experimental data, Fig. 6 illustrates a set of experi- 
mental and simulated currents measured in three ad- 
ditional photoreceptors, each stimulated at intensities 
just above threshold, near the cr value for the steady- 
state amplitude and near amplitude saturation. Tables 4 
and 5 list parametric values corresponding to simula- 
tions in Figs. 5 and 6, as well as the mean ± SD of values 
determined in a total of 10 cones. 

Increasing the intensity of light steps beyond that 
which saturates photocurrent amplitude prolongs the 
response duration, even though the stimulus duration 
is unchanged (Fig. 6 C) . To simulate step photocurrents 
beyond amplitude saturation, we adjusted lmax while 



keeping a PDf , constant at its highest 



final 



value. Fig. 6 D 



is a time-expanded view of the end of the amplitude- 
saturated step response in Fig. 6 C. Shown are simulated 
photocurrents that fit experimental data well by adjust- 
ing lmax (Fig. 6 D, black traces), as well as a simulation 
that fails to fit the experimental data (Fig. 6 D, red trace) 
because lmax was kept at the same value used to success- 
fully fit the response to dim light steps. Thus, just 
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Figure 5. Experimental (gray traces, noisy) and simulated (black traces, noiseless) photocurrents measured at —40 mV in a dark- 
adapted bass single cone. Photocurrents were elicited with 2-s light steps of intensity: 21, 61, 120, 246, 601, 1,177, 12,614, 25,873, and 
50,681 VP*/s. The values of the parameters used in the simulations shown are listed in Table 4 (Cone 1). 
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as postulated in simulations of flash photocurrents, 
the rate-limiting reaction that controls PDE* life- 
time changes from PDE* inactivation ( ) to VP* 
phosphorylation ( l max ) as a function of light is im- 
perative to successfully match simulated to experi- 
mental photocurrents. 

Under the assumptions detailed above, we success- 
fully simulated step photocurrents elicited by intensities 
between 10 and ^5 x 10 4 VP*/s, the intensity range 
tested here. The light dependence of the rates of PDE 
inactivation and VP* phosphorylation inferred from 
simulations in 10 cones is illustrated in Fig. 7 (A and B) . 
Both a. PDE and l max are light-independent at dim inten- 
sities. Above ^1,000 VP*/s, o- PDE increases exponen- 
tially from an initial value ""'a PDE up to a final final a- PDE - 
ymax is constant for all intensities below amplitude satu- 
ration; above that intensity, good fits were obtained 
keeping a. PDE at its largest final value and progressively 
reducing y max - 



Examination of Figs. 5 and 6 shows that simulations 
failed to fit the large rebound in current observed in 
the dark after the end of very bright light steps (greater 
than =25,000 VP*/s). Hamer (2000; Hamer etal., 2005) 
succeeded in simulating a similar rebound in the dark 
observed in amphibian rods by assuming there is a back 
reaction that converts excited VP* back to VP, as first 
proposed by Forti et al. (1989). It is also possible that 
the light-like response in the dark at the end of very 
bright steps arises from the ability of bleached VP to 
activate the phototransduction cascade but with much 
lower efficiency than photo-excited VP*, a phenomenon 
proposed by Matthews et al. (1996) and referred to as 
"bleaching adaptation." We do not find evidence to 
support or deny either hypothesis in bass cones, and, 
once photocurrents are saturated, simulations cannot 
be used to explore VP*, PDE*, or GC dynamics because 
changes in cGMP are no longer reflected by CNG 
channel activity. 



TABLE 4 

Values of model parameters that best simulate, photocurrents in dark-adapted bass single cones: Response to 2-s steps of light 





Parameter 


Category* 


Units 




Cone 1 






Cone 2 






Cone 3 




VP* 


Intensity 


Statistical 


VP*/s 


61 


1,177 


25,872 


171 


773 


35,236 


50 


1,975 


41,466 




Ymax 


Adjustable 


(iM/s 


100 b 


100 


90 


100 


100 


63 


100 


100 


55 


T* 


*P 

T () 


Adjustable 


1 's 


250 c 


250 c 


250 c 


250 


250 


250 


220 


220 


220 


PDE* 




Adjustable 


(iM/s per 
molecule 


0.439 d 


0.439 d 


0.439 d 


0.347 


0.347 


0.347 


0.317 


0.317 


0.317 






Adjustable 


Is 


5 e 


6.5 


2 


9 


10 




3 


5 










Is 




9 


97 






55 




8.4 


42 




b 

PDE 




Is 




9 


2 










11 




Ca 2t influx 


K m 


Adjustable 


|i.M 


0.073 


0.073 


0.073 


0.052 


0.052 


0.052 


0.098 


0.098 


0.098 




('HA 


Adjustable 


uM 


21.1 


21.1 


21.1 


28.2 


28.2 


28.2 


30.8 


30.8 


30.8 




B 


Adjustable 




11.3 


11.3 


11.3 


10 


10 


10 


15 


15 


15 


Ca 2t efflux 


Tmax 
JCa 


Statistical 


pA 


4.59 


4.59 


4.59 


5 


5 


5 


3.56 


3.56 


3.56 




exc 


Adjustable 


uM 


0.025 


0.025 


0.025 


0.025 


0.025 


0.025 


0.025 


0.025 


0.025 



'Values of invariant parameters first used to compute dark and flash photocurrents are not repeated here. They have the same values listed in Table 2. 
b J max experimental value is 100-150 s" 1 (Tachibanaki et al., 2005). 

This value was initially constrained to 150 s -1 from experimental data (Leskov et al., 2000). 

Experimental values of PDE* k cat range from 2.2 to 8 x 10 3 s -1 (Dumke et al., 1994; D 'Amours and Cote, 1999; Muradov etal., 2009). From the definition 
°f ftsub (Eq. 3.11) and the highest k cat value, the mean value (3 S1(6 0.334 {Table 5) indicates that the effective intracellular volume for cGMP is, on 
average, 0.08 pL, rather than the 0.19 pL used throughout these calculations. This factor of approximately two can probably be accounted for by the fact 
that the volume, in fact, changes some from cell to cell and is not the rigid assigned value of 0.19 pL. There is experimental uncertainty because the size 
of each and every cone recorded from was not measured. Also, the actual concentration of cGMP may be a little smaller than the computed value because 
there is intracellular cGMP binding (Olson and Pugh, 1993) , a fact that would be revealed as a smaller apparent volume. 

This value was initially constrained to 10 s -1 at dim lights from the analysis of dark current noise in bass single cones (Holcman and Korenbrot, 2005). 
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Challenging the cone phototransduction model: Predicting 
the response to flashes superimposed on light steps 
We tested the simulation model by its ability to predict 
experimental data. We investigated light adaptation, the 
observation that the response to a constant intensity 
flash superimposed on a light step becomes smaller and 
faster as the intensity of the step increases (Fig. 8). The 
flash responses are light adapted because, although 
small, they are larger than predicted from the light 
dependence of the amplitude of flash photocurrents 
measured in the same cone when it is thoroughly dark 
adapted. We first fit simulated currents to the responses 
elicited by the light steps alone, and then used the model 
to predict the response to a flash of the same intensity as 
the experimental one (Fig. 8 A) . Fig. 8 B aligns the pre- 
dicted and measured flash-elicited responses at various 
levels of background light intensity. There is good match 
between predicted and observed flash photocurrents at 
each of the various states of light adaptation tested. 

We further tested the model by challenging its ability 
to match the response to a constant, bright flash pre- 
sented at the end of step stimuli of varying intensity. Al- 
though the intensity of the bright test flash was constant 
and sufficient to saturate the photocurrent amplitude, 
the duration of the response to the flash shortened as 
the intensity of the conditioning light step increased 
(Fig. 9 A) . The paradoxical observation that the duration 
of the flash-elicited photocurrent shortens as the step 
intensity increases reproduces similar observations by 
Matthews and Sampath (2010) in tiger salamander cones. 



The simulation model successfully fit the paradoxical 
findings. Fig. 9 A presents simulated photocurrents 
optimally fit to the experimental photocurrents. Various 
experimental tests led Matthews and Sampath (2010) to 
propose that the rate-limiting reaction that controls the 
duration of the bright flash response is VP* inactivation 
(quenching). In agreement, we found that simulations 
fit experimental data only by assigning VP* phosphory- 
lation a rate-limiting role and adjusting the value of 
Imax during the flash response as a function of the in- 
tensity of the light step. Fig. 9 B illustrates the ends of 
the photocurrents activated by the same bright flash 
superimposed on the various light steps. The figure offers 
a time-expanded view to better illustrate the shortening 
of the flash-response duration as the step intensity in- 
creases. Fig. 9 C presents superimposed, time-expanded 
views of the ends of the simulated flash response that 
best fits the experimental data at each step intensity 
tested (the same data as in Fig. 9 A). The value of lmax 
that fit the step responses was the same at all intensities. 
The simulated data shown as continuous traces in Fig. 9 C 
were computed by adjusting y max for the flash response 
at each step intensity. If lmax is not adjusted but simply 
kept at the value used to compute the step response, 
simulations do not fit the experimental data (for exam- 
ple, dashed line in Fig. 8 C). The absolute need to ad- 
just the lmax value to simulate the experimental results 
affirms that the dominant rate constant that controls 
response duration at intensities that saturate photocur- 
rent amplitude and up to ^5 x 10 4 VP*/s (the range 



TABLE 5 

Values of model parameters that best simulate photocurrents in dark-adapted bass single cones: Response to 2-s steps of light 





Parameter 


Category* 


Units 




Cone 4 




Mean ± SD 


VP* 


VP*/s 


Statistical 


VP*/s 


50 


1,975 


41,466 


Fig. 7 




"lmax 


Adjustable 


uM/s 


95 


95 


40 




T* 




Adjustable 


1/s 


250 


250 


250 


237 ± 31 


PDE* 




Adjustable 


pM/s per 
molecule 


0.239 


0.239 


0.239 


e 0.366± 0.101 






Adjustable 


1/s 


2.7 


2.8 




Fig. 7 




final 

PDE 




1/s 




7.8 


55 






h 

PDE 




1 s 




10 






Ca 2+ influx 


K HA 


Adjustable 


uM 


0.12 


0.12 


0.12 


0.085 ± 0.024 




Cm 


Adjustable 


uM 


31.4 


31.4 


31.4 


26.4 ± 9.3 




B 


Adjustable 




15 


15 


15 


11 ± 3.3 


Ca 2+ efflux 


Tmax 
JCa 


Statistical 


pA 


3.89 


3.89 


3.89 


3.96 ± 0.9 




exc 
Kca 


Adjustable 


uM 


0.030 


0.030 


0.030 


0.025 ± 0.006 



Mean ± SD; n - 10 cones. 

a Values of invariant parameters first used to compute dark and flash photocurrents are not repeated here. They have the same values listed in Table 2. 
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explored here) is the Ca-dependent rate of VP* inacti- 
vation by phosphorylation. 

Computational gene knockout: The functional role of cone 
CNG ion channel modulation 

The molecular model of cone phototransduction de- 
tailed here incorporates Ca-dependent modulation of 
the cGMP sensitivity of ion channels that sustain the 
photocurrent. The quantitative features of this modula- 
tion are known in detail in bass single cones (Rebrik 
and Korenbrot, 1998; Rebrik et al., 2000), as is the iden- 
tity of a novel soluble Ca 2+ -binding protein named CNG- 
modulin that mediates this effect (unpublished data) . 
The specific physiological role of channel modulation 
in the phototransduction process, however, is unknown. 
We investigated what this role might possibly be by 
examining with photocurrent simulations the conse- 
quence of deleting channel modulation, what might 
be called "computational gene knockout." 

We will refer to the model detailed above as "com- 
plete" and as "CNG channel modulation-minus" the 
model in all equations identical to the complete one, ex- 
cept Ca-dependent channel modulation is omitted by re- 
placing Eqs. 2.5 and 2.6 by a single equation (Eq. 2.5) with 
unchanging KcGMP =172 uM, the value correspondent 



to Ca 2+ = 0.4 pM. Figs. 10 and 11 show side-by-side com- 
parisons of current simulations computed with the 
complete and the CNG channel modulation-minus 
models in response to flash or step stimuli of varying 
intensity. In each figure, the left panels illustrate simula- 
tions computed with the complete model, and the right 
panels illustrate those computed with the CNG channel 
modulation-minus one. Similar effects were observed in 
all other sets of experiment/ computational data tested. 

The flash photocurrents and simulated data shown in 
Fig. 10 are a subset of the data measured in cone 1 of 
Fig. 2. In the absence of CNG channel modulation, the 
initial rate of rise of flash photocurrents is hardly 
affected because VP* and PDE* activation are not Ca 
dependent. But photocurrent peak sensitivity is en- 
hanced and its time integral increased, and the "off 
response is much larger in amplitude and the first one 
in a sequence of dampened oscillations at the end of 
the response at all light intensities. 

Increased photosensitivity and oscillations in the 
photocurrent are experimentally observed when Ca 2+ - 
buffering molecules (BAPTA, for example) are loaded 
into the cone cytoplasm (Matthews etal., 1990; Holcman 
and Korenbrot, 2005) , but the effects of deleting CNG 
modulation are not simply the same as those of 
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Figure 6. Experimental (gray traces, noisy) and simulated (black traces, noiseless) step photocurrents measured at —40 mV in three 
different dark-adapted bass single cones. (A) Intensities tested were 171, 773, and 8,770 VP*/s. (B) Intensities tested were 50, 963, and 
10,320 VP*/s. (C) Intensities tested were 50, 963, 10,320, and 41,466 VP*/s. The values of the parameters used in the simulations shown 
are listed in Table 4. (D) A time-expanded view of the end of the response to 41,466 VP*/ s shown in C. At this intensity, photocurrent 
amplitude was saturated. The unsuccessful simulations illustrated in red were computed when Ymax was kept at the same value used to 
fit the responses to dim stimuli. 
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"enhanced" buffering: Buffering prolongs response 
duration to a much larger extent than does deleting 
channel modulation, and the oscillations at "off' are 
much smaller than those anticipated by neglecting 
CNG channel modulation. It stands to reason that at- 
tenuating light-dependent changes in cytoplasmic Ca 2+ 
by BAPTA would have overlapping effects with those an- 
ticipated when channel feedback is deleted; buffers 
attenuate by reducing the magnitude of Ca 2+ concentra- 
tion changes, whereas deleting the channel attenuates 
by removing a target of Ca 2+ effect. Moreover, Ca 2+ buf- 
fers affect more than channel modulation; they also 
affect the other Ca-dependent reactions in the photo- 
transduction cascade (VP* phosphorylation and GC 
activity). Hence, buffers can affect functional features 
not affected when channel modulation is deleted. 

The step photocurrent and simulated data shown in 
Fig. 11 are a subset of data measured in yet another 
cone, not illustrated in prior figures. In the absence of 
CNG modulation, the photocurrent initial rate of rise is 
hardly changed. More remarkable, the step response 
becomes much less stable. At dim light steps, promi- 
nent dampened oscillations are anticipated both dur- 
ing and after the light step. At brighter light steps, 
photocurrents are more light sensitive at both their 
peak and steady state, and light adaptation is compro- 
mised, as gauged by the lessened sag from peak to 
steady state. In summary, Ca-dependent CNG channel 



modulation in cones introduces a critical regulatory 
function, a quality control function, which assures that 
the complex sequence of chemical reactions underly- 
ing phototransduction yield the desired speed, light 
sensitivity, and light-adaptation features, and operate with 
great stability. 

DISCUSSION 

Mathematical models offer a succinct and precise ex- 
pression of the state of understanding of a given physi- 
ological process at the time the model is developed. We 
present here new experimental data along with an 
explanatory model of cone phototransduction that ex- 
pands prior cone transduction models to include explicit 
characterization of Ca-dependent VP phosphorylation 
and Ca-dependent CNG channel sensitivity modulation. 
Using the model, we simulate bass cone photocurrents 
elicited by light flashes, steps, and flashes superimposed 
on steps and fit the experimental results well. 

We measured photocurrents in isolated bass single 
cones using tight-seal electrodes in the whole cell mode. 
Two technical details merit highlighting. First, this is 
the first study we know of to measure photocurrents 
elicited by light flashes, steps, and flashes superimposed 
on steps in the same cone under voltage clamp using 
tight-seal electrodes. This is important because in the 
absence of voltage clamp, the temporal dynamics of 
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Figure 7. Light dependence of the values of ad- 
justable parameters inferred from successful simu- 
lations of step photocurrents. Data points are the 
mean ± SD of values computed in 10 single cones. 
The light intensity shown on the graph is the mean 
of intensities binned within plus or minus 10% of 
the mean. (A) Initial and final values of PDE* inacti- 
vationrate ( intt a Pm and ""a fDE ; Eq. 3.14). (B) VP* 
phosphorylation rate ( Imax ) . The continuous lines 
are drawn by eye to join the data points. 
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photocurrents measured in large cone outer segments 
are distorted because of the large membrane capaci- 
tance (in the range between 65 and 100 pF) (Miller and 
Korenbrot, 1993) . Second, there is a serious disadvan- 
tage in the use of tight-seal electrodes in the whole cell 
mode because small-sized, soluble cytoplasmic mole- 
cules are lost by diffusion into the electrode lumen, and 
the cytoplasmic ionic composition may be different 
than in the intact cell, including the features of the 
cytoplasmic Ca 2+ buffer. Indeed, in previous studies of 
bass single cones, we have noted the change in photo- 
current kinetics upon attaining whole cell mode (Miller 
and Korenbrot, 1993) and the slow loss of Ca-dependent 
CNG ion channel modulation, which we attributed to 
the diffusion loss of the modulator molecule (Rebrik 
and Korenbrot, 1998, 2004). We mitigated this prob- 
lem by limiting data reported here to photocurrents 
measured within the first 7-8 min after attaining 
whole cell mode, an interval over which photoresponse 



kinetics are stable and reproducible, although they can 
change thereafter. 

The model makes explicit assumptions about the 
features of intracellular Ca 2+ buffers in cone outer seg- 
ments, about which very little is actually known. Two 
different models have been used previously to simulate 
the buffer in photoreceptors. The first one, most fre- 
quently used, assumes that buffer molecules bind Ca 2+ 
instantaneously in the time scale of events in photo- 
transduction, and the buffer is simply characterized by 
its capacity, the ratio of bound/free Ca 2+ (Miller and 
Korenbrot, 1994; Soo et al., 2008). We found this as- 
sumption inadequate in the simulations reported here 
because it does not allow for time-dependent changes 
in buffer capacity. The second model was introduced by 
Forti et al. (1989) and assumes that buffering capacity is 
time dependent because Ca 2+ -binding molecules become 
available in a time-dependent manner. This model was 
developed to explain the sag in the step response of rod 
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Figure 8. Experimental (gray traces, noisy) and simulated (black traces, noiseless) photocurrents elicited by the same intensity flash 
superimposed on light steps of varying intensity. (A) Photocurrents generated by a 150-VP* flash delivered 1 s after the onset of light 
steps of intensities: 171, 773, 8,770, and 35,236 VP*/s. (B) Time-expanded view of the response to the constant intensity flash. Flash 
photocurrent on 171 VP*/ s background was 7.3 pA in peak amplitude, 80 msec time to peak. On 773 VP*/s background it was 6.5 pA 
in peak amplitude, 77 msec time to peak. On 8,770 VP*/ s background it was 2.9 pA in peak amplitude, 73 msec time to peak. On 35,236 
VP*/s background it was 0.9 pA in peak amplitude, 63 msec time to peak. 
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photoreceptors and has been extensively applied by 
Hamer and coworkers in their studies of rod phototrans- 
duction (Hamer, 2000; Hamer et al., 2005). Under this 
assumption, optimal fit of simulated data to experimental 
data in rods predicts that light causes rapid decreases in 
cytoplasmic free Ca 2+ , 50% of which are completed within 
tens of milliseconds from light onset. But this predicted 
fast component is not observed in experimental measure- 
ments in rods (Ratto et al., 1988; Gray-Keller and Detwiler, 
1994) or cones (Sampath et al., 1999; Leung et al., 2007) . 
For this reason, we elected not to adopt this assumption. 

The model we adopted allows time dependency be- 
cause the buffer capacity is not constant but a function 
of free Ca 2+ ; buffer capacity changes in time as Ca 2+ con- 
centration changes. The quantitative features of a high 
affinity buffer inferred from simulations are close to 
values experimentally measured in other cells, includ- 
ing rods (Lagnado et al., 1992; Neher and Augustine, 



1992; Berlin et al., 1994; Schwiening and Thomas, 1996; 
Xu et al., 1997; Trafford et al., 1999). The adequacy of 
the values of Ca 2+ -buffering parameters arrived at by 
simulation is best justified by the fact that the rate of 
Ca 2+ clearance from bass cone outer segments, the time 
to reduce Ca 2+ by 1/e when all CNG channels are closed 
by light, is, on average, 40 msec in simulations, essen- 
tially the same value measured experimentally (Sampath 
etal., 1999). 

The duration of the photoresponse in cones is deter- 
mined by the lifetime of light-excited PDE*, a feature 
dominated by the rate constant of one of three molecu- 
lar events, depending on stimulus intensity: PDE* inac- 
tivation ( apDE ) , VP* inactivation by Ca-dependent 
phosphorylation ( Xmax ), and the rate of decay of the 
cone opsin photoproduct Metall. The only hypothesis 
that allowed for successful simulation of all patterns 
of stimulation (flash, step, and flash superimposed on 
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Figure 9. Experimental (gray traces, 
noisy) and simulated (black traces, 
noiseless) photocurrents elicited by 
the same intensity flash superimposed 
on light steps of varying intensity. 
(A) Photocurrents generated by the 
same flash, one that saturated photo- 
current amplitude, delivered at the end 
of 2.5-s light steps of intensities: 171, 
773, 8,770, 35,236, and 206,056 VP*/s. 
(B and C) Time-expanded views of 
experimental and simulated photo- 
currents, respectively. Shown are the 
ends of the saturated flash responses 
recorded at various step intensities, as 
labeled. Successful flash simulations 
were obtained by adjusting the value 
of ^ max from 90 s _1 during the step 
response to 90 s _1 at 171 VP*/s, 90 s -1 
at 773 VP*/s, 52 s" 1 at 8,770 VP*/s, 
28 s -1 at 33,265 VP*/s, and 30 s _1 at 
206,056 VP*/s. The dashed black line 
in C is the simulated response on a 
background of 206,056 VP*/ s com- 
puted without adjusting ~^ max ' keep- 
ing its value the same as that used to 
simulate the step response. 
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step) and all intensities postulates that the dominant 
rate constant that controls PDE* lifetime changes with 
light intensity: it is Upde at intensities below amplitude 
saturation (less than ^4,000 VP* in a flash; M).001% VP 
bleached) and it is l ma x at intensities above that and up 
to ^5 x 10 4 VP* (M).01% VP bleached). When intensi- 
ties bleach 0.1% of the VP and above, the rate of Mil decay 
is the dominant rate constant (Estevez et al., 2009). 

The postulated light-dependent change in the domi- 
nant rate constant of photocurrent duration is consis- 
tent with findings in cone photoreceptors and is in 
contrast to the understanding of the dominant rate 
constant in rods (Burns and Pugh, 2009). For flashes 
activating up to ^4,000 VP* in rods, the dominant 
rate constant is « PD£ , the rate of association of the RGS9 
multiprotein complex with the T*-PDE* complex. 



For flashes activating 4,000-160,000 VP*, the dominant 
rate constant is the rate of association of RGS9 with T* 
alone, molecules activated by VP* but not forming a 
complex with PDE* (Burns and Pugh, 2010). In cones, 
and over this range intensity, in contrast, the dominant 
rate constant controlling photocurrent duration is Ymax, 
the rate of VP* phosphorylation. 

A feature in our model, not included in any previous 
models, is the Ca-dependent modulation of the cGMP 
dependence of CNG channel gating. The most signifi- 
cant consequence of this modulation is that in cones, 
CNG ion channel activity is controlled by both cGMP 
and Ca 2+ . cGMP controls the probability of channel 
opening (gating) (Picones and Korenbrot, 1994), whereas 
Ca 2+ controls the effectiveness of a given cGMP concen- 
tration to gate (Rebrik etal., 2000). Our computational 
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Figure 10. The physiological role of Ca- 
dependent CNG channel modulation in 
flash photocurrents. Experimental photo- 
currents were measured at —40 mV in a 
dark-adapted bass single cone in response 
to 10-msec flashes of intensities: 71, 167, 
710, and 3,561 VP*. The panels on the left 
illustrate experimental (gray traces, noisy) 
and simulated (black traces, noiseless) 
photocurrents computed with the normal 
complete model (Table 3, Cone 1). The 
panels on the right reproduce the same 
experimental data (gray traces, noisy) , now 
superimposed by simulated photocurrents 
(black traces, noiseless) computed with a 
model in which CNG channel modulation 
is omitted. 



50 Modeling photocurrents in cones 



simulations reveal that no one single specific feature of 
the phototransduction signal is controlled by the chan- 
nel modulation. Rather, just as Ca 2+ itself, ion channel 
modulation introduces one more regulatory tool to 
achieve the speed, light sensitivity, adaptation, and sta- 
bility characteristic of the cone photoresponse. 

CNG channel modulation is a not a prerequisite for 
stability in the phototransduction enzymatic cascade. It 
appears necessary only for the relatively fast dynamics of 
cone phototransduction. Rod photoresponses are stable, 
even though in these cells CNG channel modulation is 
very small in extent and occurs only at very low cytoplas- 
mic Ca 2+ levels (tens nanomolar) (Nakatani et al., 1995; 
Bauer, 1996; Sagoo and Lagnado, 1996). In fact, knock- 
ing out channel modulation in mammalian rods is of lit- 
tle functional consequence (Chen etal., 2010). 



To understand why CNG ion channel modulation is 
of such functional importance, we illustrate with simula- 
tions some of the molecular events associated with flash 
and step responses in the complete and the CNG 
channel modulation-minus models. In both flash and 
step responses, the initial rate of change of the outer 
segment inward currents in flash and step photore- 
sponses is essentially the same in the presence or ab- 
sence of ion channel modulation (Fig. 12 A). This is 
because the activation of the VP* and PDE*, the feature 
that determines the initial rate of photocurrent rise, is 
not Ca dependent. 

Compared with the complete model, absence of 
channel modulation causes a larger change in current 
amplitude and a prolonged time integral in response to 
flash stimulation (Fig. 12 A) . Because of these differences 
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Figure 1 1 . The physiological role of Ca- 
dependent CNG channel modulation in 
step photocurrents. Experimental photo- 
currents were measured in a dark-adapted 
bass single cone in response to steps of 
61, 124, 1,177, and 2,414 VP*. The panels 
on the left illustrate experimental (gray 
traces, noisy) and simulated (black traces, 
noiseless) photocurrents computed with 
the normal complete model (Table 4, 
Cone 1). The panels on the right repro- 
duce the same experimental data (gray 
traces, noisy), now superimposed by sim- 
ulated photocurrents (black, noiseless 
traces) computed with a model in which 
CNG channel modulation is omitted. 
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in the photocurrent, there is a difference in the light- 
dependent changes in cytoplasmic Ca 2+ . In the com- 
plete model, Ca 2+ decreases and smoothly returns to 
its starting level, the same as is experimentally ob- 
served (Leung et al., 2007). In the absence of chan- 
nel modulation, cytoplasmic Ca 2+ falls to a lower level 
than normal and stays at that lower level for a longer 
period of time (Fig. 12 B). These seemingly small dif- 
ferences have enormous consequence because GC 
activity is so sensitive to free cytoplasmic Ca 2+ : GC ac- 
tivity is enhanced to a greater extent and over a lon- 
ger time period in the absence of channel modulation 
than in its presence (compare GC activity in Fig. 12, 
C and D). 

In the model of a normal cell, the dynamic balance be- 
tween the Ca-enhanced GC activity and the light-enhanced 



PDE activity results in a properly controlled decrease 
and recovery of cGMP concentration. PDE activity ini- 
tially exceeds that of GC and then lags behind it until 
the activity of both enzymes returns to the same initial 
dark value (Fig. 12 C, flash response). In the absence 
of channel modulation, this dynamic balance is lost 
(Fig. 12 D). PDE activity is essentially unchanged, but 
the enhanced GC activity causes a larger increase in 
cGMP concentration and, therefore, larger changes in the 
inward membrane current in the CNG modulation-minus 
model than the complete one (Fig. 12 A). This en- 
hanced inward current results in cytoplasmic Ca 2+ over- 
load (Fig. 12 B), which, in turn, completely shuts down 
GC activity (D) . Although the PDE activity is small be- 
cause it has returned to its dark level, it is nonetheless 
higher than that of GC (Fig. 12 D). Hence, cGMP 
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Figure 12. Simulated dynamics of 
several molecular events underlying 
flash and step photoresponses in bass 
single cones. Left panels show various 
biochemical and biophysical events 
elicited by a 10-msec flash of 167 VP* 
intensity (Table 3, Cone 1). Right pan- 
els show the same events elicited by 
a 2-s step of 124 VP*/s in a different 
cone (Table 4, Cone 1 ) . (A and B) The 
inward outer segment membrane cur- 
rent and the cytoplasmic free Ca 2+ 
concentration. Superimposed in each 
panel are the results of simulations 
with the complete model (blue traces) 
and one in which CNG channel modu- 
lation is omitted (red traces). (C) En- 
zymatic activity of PDE (dashed line) 
and GC (solid line) and their change 
with light computed with the complete 
model. (D) Enzymatic activity of PDE 
(dashed line) and GC (solid line) com- 
puted with the model that omits CNG 
channel modulations. 
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concentration decreases, channels close, and the anom- 
alously enhanced inward current returns toward its 
starting value. As the inward current decreases, so does 
the cytoplasmic free Ca 2+ (Fig. 12 B). The dynamic of the 
PDE-GC balance, however, is not at equilibrium; the 
iteration of over/under cGMP synthesis is repeated until 
a true equilibrium is reestablished, thus the oscillations 
in the membrane current. 

In the complete model, the response to a steady light 
attains a stationary value because the activity of PDE and 
GC reach a new balance and are the same as long as light 
is present (Fig. 12 C). Again, the most dramatic conse- 
quence of the lack of channel modulation is a change in 
the extent and time course of Ca 2+ concentration changes, 
which, in turn, change the time course and extent of 
changes in GC activity. Because PDE activity is little 
changed by the lack of channel modulation, the normal 
dynamic balance of the activity of the two enzymes is 
severely disturbed (Fig. 12, compare C and D, right 
column) . The enzymatic activities are not in balance, and 
cycles of over/ under cGMP synthesis are repeated, caus- 
ing oscillation in membrane current both while the light 
is present and at the termination of the light step. 

In brief, the simultaneous control of CNG ion channel 
activity by cGMP and Ca 2+ causes the kinetics of light- 
dependent channel closing and reopening to be far dif- 
ferent than they would be if controlled by cGMP alone. 
Under this dual control, the light-dependent changes in 
cytoplasmic free Ca 2+ and, therefore, changes in GC ac- 
tivity are of appropriate amplitude and time course and 
in dynamic balance with PDE activity. If channels were 
controlled by cGMP alone, the light-dependent changes 
of Ca 2+ would be far different, and the balance between 
light-dependent changes in GC and PDE activities would 
be lost. 

Ca-dependent control of CNG channel ligand sensitiv- 
ity and of VPK activity is an experimental fact that must 
be included in complete models of cone phototransduc- 
tion. Yet, a contemporary model of cone phototransduc- 
tion (Soo et al., 2008) that ignores Ca-dependent control 
of VP* phosphorylation and channel modulation can 
simulate photocurrents that well fit experimental data. 
In our own exploration, we find that even if the Ca- 
dependent regulation is ignored, judicious selection of 
values for the various adjustable parameters could 
simulate some experimental data, although for some 
parameters, their values are often different than the ex- 
perimentally known value, and the models do not work 
nearly as well over the many stimulation paradigms and 
intensities reported here. This only indicates that what is 
gained by Ca-dependent controls is not the difference 
between transduction and the failure to transduce, but 
the ability to regulate the molecular events underlying 
phototransduction to attain the biological imperatives of 
time course, sensitivity, and light adaptation that sustain 
the visual process. 
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